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Abstract 
Properties of nuclear and neutron matter are discussed in a nonlinear a-to-p mean-field approxima- 
tion with self- interactions and mixing-interactions of mesons and baryons. The nonlinear interactions 
are renormalized by employing the theory of conserving approximations, which results in a thermo- 
dynamically consistent approximation that maintains Hugenholtz-Van Hove theorem and Landau's 
requirement of quasiparticles. The approximation is equivalent to the Hartree approximation with 
effective masses and effective coupling constants of baryons and mesons. The effective masses and cou- 
pling constants are naturally required by self-consistency of the theory of conserving approximations. 
The approximation is applied to nuclear and neutron matter, which suggests that the lower bound of 
nuclear compressibility K ~ 180 MeV (with the symmetry energy 04 = 35.0 MeV) be required to be 



consistent with properties of nuclear matter and the maximum masses of observed hadronic neutron 
stars (M max > 2.00 Mq). The values of the compressibility, symmetry energy together with effective 
masses and coupling constants of baryons and mesons will be important constraints to examine mod- 
els of nuclear and neutron matter. The accumulating data and accurate measurements of observables 
in high density and energy region will supply significant information in order to testify theoretical 
consistency of nuclear models. 
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1 Introduction 

The linear a-to mean-field approximation has been extensively applied to finite nuclei, nucleon-nucleon 
interactions as well as to equations of state of nuclear matter and properties of neutron stars ' ' >. 
The structure of the linear mean-field approximation is extended to nonlinear mean-field models which 
have also been extensively applied to properties of nuclear and neutron matter!' 1 '' '. Neutron stars 
have attracted much interest in nuclear and high energy physics with a large body of accumulating 
experimental evidence. The observed masses of hadronic neutron stars are above 1.3 Mq (the solar 
mass: M ~ 1.989 x 10 30 kg), and the maximum masses of neutron stars are expected to be below 
2.5 MqI ', or 2.2 ± 0.2Mq' >. Since a neutron star is mostly composed of hadrons, the properties 
of a neutron star may be calculated with the self-consistent equation of state for the many-body 
system of hadrons. We will employ a nonlinear mean-field model to investigate relations between 
properties of nuclear matter and neutron stars constrained by the maximum masses of neutron stars: 
M max = 2.50 M Q , and 2.00 M . 

Approximations for infinite nuclear and neutron matter are required to maintain self-consistent 
conditions: the Hugenholtz-Van Hove theorem *■ ' (HV theorem), the virial theorem' > and the theory 
of conserving approximations' >. The self-consistent conditions are essential to calculate explicitly 
microscopic and macroscopic physical quantities and denoted as thermodynamic consistency' ' ' 
to an approximation. Serot and Walecka's approach is extended by maintaining thermodynamic 
consistency in the current nonlinear model, and the extension leads to the concept of effective masses 
of mesons and effective coupling constants. The mean-field approximation with nonlinear interactions 
of baryons and mesons will become a conserving approximation when effective masses of mesons 
and effective vertex coupling constants are self-consistently defined. This results in an important 
conclusion that the nonlinear interactions of mesons should be examined as density-dependent, or 
energy-dependent effective masses and coupling constants. 

The nonlinear mean-field model contains several adjustable parameters coming from self-interactions 
and mixing- interactions of mesons and baryons. Typically, the treatment of quantum vacuum is ne- 
glected entirely that it may be difficult to study nuclear systems at short distances, or high-energy 
regions' ' '. However, the current nonlinear mean-field approximation constructed with thermody- 
namic consistency is a self-consistent approximation which constrains nonlinear interactions and ad- 
justable parameters in terms of physical conservation laws. In this sense, the general features obtained 
in the calculation may compensate for the deficiency and remain valid as model-independent results. 

The nonlinear mean-field approximation will be proved equivalent to Hartree approximation with 
effective masses and effective vertex coupling constants. The results indicate that nonlinear interac- 
tions have a significant effect on properties of nuclear and neutron matter; however, since nonlinear 



interactions are strongly confined by self-consistency, their contributions to physical quantities are 
obtained in a balanced and restricted way. Hence, the existence and strength of nonlinear interactions 
of hadrons could be examined if effective masses of mesons are measured from the nuclear experi- 
mental data. The renormalized effective masses and vertex coupling constants are derived explicitly, 
and physical contributions of nonlinear interactions to properties of nuclear and neutron matter are 
checked quantitatively. The analysis of nonlinear interactions by way of effective masses and coupling 
constants is one of the purposes of the current paper. The self-consistency of meson-propagators with 
effective masses of mesons and effective coupling constants will be explained in detail in the sec. 3 and 
sec. 4. 

In order to simulate complicated interactions and retardation effects, the nonlinear a mean-field' ' 
and the nonlinear a-to-p mean-field models' ' have been discussed. By adjusting coupling constants 
of nonlinear self-interactions of mesons, it seems that one could obtain softer equations of state as 
one wishes' > in order to simulate nuclear properties. However, since the nonlinear interactions can 
be renormalized as effective masses and coupling constants, contributions to physical quantities will 
be restricted by self-consistency; hence, there appear a possibility to check the validity of nonlinear 
interactions by experimental data of nuclear matter and neurton stars. Furthermore, if nonlinear 
mixing interactions are introduced, coupled equations of motion for mesons and baryons will confine 
the values of coupling constants more strictly. Therefore, coupling constants of nonlinear interactions 
are not free parameters at all. Coupled equations of motion and self-consistency as a conserving 
approximation will restrict the strength of nonlinear interactions. This is an important consequence 
derived from the current conserving mean-field approximation. 

The properties of neutron stars, such as the mass, moment of inertia and radius of neutron stars, 
are calculated by employing the equation of state in the nonlinear mean-field approximation. The 
coupling constants are confined within certain values in order to produce saturation properties and 
the symmetry energy, a^ = 35.0 MeV' >, of nuclear matter and masses of neutron stars. We have 
calculated properties of nuclear matter and neutron stars (nuclear compressibility, symmetry energy, 
effective masses of nucleons and mesons, the maximum mass of neutron stars) and investigated the 
minimum value of nuclear compressibility, K, to the value of the supposed maximum mass of neutron 
stars, M max = 2.50 M & and 2.00 M Q . The nonlinear mean-field approximation suggests that a lower 
bound of compressibility K > 340 MeV should be required to be compatible with the maximum value, 
M star > 2.50 M ; and K > 180 MeV compatible with M star > 2.00 M . 

It is difficult to reproduce K and 04, M max , simultaneously by adjusting the values of coupling con- 
stants. The conditions of thermodynamic consistency and coupled equations of motion for mesons and 
baryons strictly restrict contributions of nonlinear interactions. Consequently, the values of nonlinear 



coupling constants are confined in a narrow region. Therefore, nonlinear interactions will generate 
certain fixed experimental values, and the validity of nonlinear interactions and relations between 
properties of nuclear matter and neutron stars may be examined by comparing with experimental 
data of nuclear matter and neutron stars. These results are explained in the sec. 5. Conclusions and 
remarks on the self-consistent nonlinear mean-field approximation are in the sec. 6. 

2 The nonlinear a-cu-p mean-field lagrangian with mixing and self- 
interactions of mesons and baryons 

The mean-field lagrangian of the nonlinear a-cu-p self-interactions, mixing- interactions, atoN and apN 
couplings of mesons and nucleons is defined by, 
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where the meson-field operators are replaced by expectation values in the ground state: cfto for the 
cj-field, Vo for the vector-isoscalar w-meson, V^V^ = Vq — V 2 , (/i = 0,1,2,3); the neutral /9-meson 
mean-field, Ro, is chosen for T3-direction in isospin space. The density-dependent effective vertex 
interactions of g a ujN an d g a pN terms are defined by way of the scalar mean- field, 4>q. Let us notice that 
the replacement of the meson quantum fields with classical fields means that the nonlinear mean-field 
lagrangian should be understood as an effective model to simulate complex many-body interactions 
as simple as possible in the mean-field approximation. 

The nonlinear model is motivated by preserving the structure of Serot and Walecka's linear a-uo 
mean-field approximation! ' as explicit as possible, Lorentz-invariance and renormalizability, thermo- 
dynamic consistency: the Hugenholtz-Van Hove theorem!"!, the virial theorem^, and conditions of 
conserving approximations [' '. The concept of effective masses and effective coupling constants are 
naturally generated by self-consistent conditions and will be proved by the requirements of theorems 
and conditions of conserving approximations. 

The masses in (2.1) are: M = 939 MeV, m a = 550 MeV, m^ = 783 MeV and m p = 770 MeV, 
in order to compare the effects of nonlinear interactions with those of the linear a-uj approximation 
discussed by Serot and WalecW 1 . The coupling constants may be supposed to be free parameters at 
the outset, however, they are constrained by saturation properties: the binding energy —15.75 MeV 
at k F = 1.30 fm _1 with the symmetry energy, 04 = 35.0 MeV, and in addition, by investigating the 
lower bound of nuclear compressibility, K, which corresponds to maximum masses of neutron stars. 
The maximum masses of observed neutron stars are chosen in the current analysis as M max = 2.50 
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M Q and 2.00 M l 3 ' "J, in order to examine contributions of self-consistently constrained nonlinear 
interactions in the level of Hartree approximation. 
The equation of motion for baryons is given by, 

iil^ - gllvVv - ^-lonRo) - (M - g^ )] if = , (2.2) 

where g* = g w + g acJ N4>o and 9* p = 9p + ^9upN<t>o are effective coupling constants for oj and p mesons, 
respectively. The equations of motion for mesons in the mean-field approach are given by, 

™^o + 9 -fV$ + ^^Fo + 9 -fR 2 V = g^H) = glp B , 

m 2 p Ro + ^«o 3 + 9 -f^Ro + f^o = f <VW) = f Ps • (2-3) 

The scalar density, p s , is defined by 

(2vr)3 7 £*(fc) 

where M* = M — g a <pQ. The baryon density will be denoted as p B = (ip^ip) = — tt (C = 4 for nuclear 
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matter; (" = 2 for neutron matter), and p3 = (ip'r^ip) = (kp — kp )/37r 2 ; the Fermi-momentum kF p is 
for proton and kp n for neutron. In the linear a-u; mean-field approximation, the scalar source is only 
generated by the scalar density, p s . However, one should note that the nonlinear mixing interactions, 
such as ^7oVoV' and if^foT^Roif j generate additional source terms to the equation of motion for the 
scalar field, <pQ. The equations of motion for uj and p mesons are formally obtained by replacing g w 
and g p in the linear mean-field approximation to g* and g*. 

The theory of conserving approximations requires that the equations of motion be derived as the 
functional derivative of energy density with respect to self-energies 11 (or equivalently, with respect 
to mean-fields: <fio, Vq, Ro), and moreover, it requires that the equations of motion expressed in terms 
of self-energies be derived by Feynman diagrams by applying propagators of baryons and mesons. 
It is essential for the theory of infinite matter to prove the equivalence of both calculations. The 
propagators of baryons and mesons are defined in the sec. 3; the functional derivative of energy 
density with respect to mean-fields of mesons and the equivalence of both approach will be proved in 
the sec. 4. 

3 The nonlinear a-ui-p Hartree approximation 

The approximation-scheme of Serot and Walecka's linear a-u approximation is maintained by em- 
ploying effective coupling constants and effective masses, and it is readily proved that the current 



approximation is a self-consistent, relativistic Hartree approximation as discussed in the linear a-oo 
approximation . 

The Green's function for baryons satisfies Dyson- Schwinger equation, and the formal solution to 
the Green's function is: 

[G(k)}- 1 = 7m (F + E"(fc)) - (M + E»(fc)) , 

= l0 (k° + £°(fc)) - 7 • k(l + E w (fc)) - (M + E s (fc)) , (3.1) 

where the spatial part of the self-energy, T, l (k), is expressed as £*(&;) = k l T, v (k) (i = 1,2,3); however, 
in the level of the mean-field approximation, the self-energy of the spatial part will vanish because of 
the spatial homogeneity, resulting in T,'"(k) = 0. As explained in the sec. 1, the quantum vacuum is 
entirely neglected, and the Green's function for the Fermi-sea particles in the rest frame of nuclear 
matter is used: 

G D (k) = ( 7 "fe* + M*(k))^rS(k° - E{k))d{k F - |k|) , (3.2) 

where k*° = k° + S°. The effective nucleon mass is M*(k F ) = M — g a 4>o = M + S s (/c F ), and the 
baryon single particle energy is E(k) = E*(k) - £°(fc), with E*(k) = (k 2 + M* 2 ) 1/2 . 

The mean-fields of mesons and self-energies of baryons are specifically associated by way of the 
formal solution to Green's function (3.1) and the equation of motion for baryons (2.2); the w-meson 
contribution to the baryon self-energy is denoted by EjJ, = — 5C^b> an d />meson contributions to the 
self-energies for proton and neutron are denoted as T,° pp = —(g*/2)Ro and E° n = (g*/2)Ro, respectively. 

According to the linear Hartree approximation! >, we will define the free, a, to and p, meson- 
propagators, D^\k), D^^ u (k), and D %b{k), for the nonlinear mean-field approximation as: 

a2°(*) = T2 W-' D^(k)= ~ 9 " U , (3.3) 

k^ — my- + is ki — mlf + is 

and 

where to*, to* and to* are effective masses of a, uj, p mesons, respectively. The effective masses are 
density-dependent constants, and therefore, they are determined when energy density, pressure and 
self-energies are self-consistently determined. 

The self-energies are calculated by applying G D (k), D^(k), Dj° ) > J ' u (k) and D^^ ' (k) upon (tadpole) 
Feynman diagrams as shown in Fig. la and Fig. lb. The interaction-lines in Fig. la and Fig. lb are 
given by effective masses of mesons which should be self-consistently defined. The vertex factor of 
baryon and p-meson coupling is —i(g*/2)T3 / y^. Let us notice that the scalar self-energy, S s , is generated 
by three density sources in the nonlinear model. The scalar self-energy given by the Feynman diagram 



in Fig. la is 
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where p* is the modified scalar density: p* = (p s — — Vop B ^— RqPz)- The w-meson and p-meson 

9a 9a 

contributions to the self-energy are given by 
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and 

E P<5) = **§» /^(^»M) = *4§^"° • < 3 - 7 > 

where T3 is the isospin matrix. Self-energies, (3.5) ~ (3.7), are formally identical with the equations 
in the linear mean-field approximation! '. As the effective mass of baryons is defined so that Dirac 
equation of motion is preserved, the effective masses of mesons are defined so that they are formally 
identical with the equations of motion of mesons in the linear mean-field approximation. 

The self-energies of mesons, (3.5) ~ (3.7), derived from Feynman diagrams are consistent with 
equations of motion of mesons (2.3). Since the self-energies and mean-fields are related with T, s (k F ) = 
-g a <Po, S °(M = -giyo, and H°,pJk F ) = =f(#*/2)_R ((-) for proton and (+) for neutron), the self- 
energies from (3.5) to (3.7) must be the same as the meson equations of motion, (2.3), respectively. 
Therefore, by comparing the self-energies from (3.5) to (3.7) and equations of motion for mesons (2.3), 
one can examine that the current nonlinear mean-field approximation will be self-consistent, only if 
effective masses, m* , m* and m* , are given by, 

,*2 „2 f-t i 9a3 , 9aA ,2 9au> ^ r 2 , C ~ \ 9ap p 2 
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m » ~ M 1 ~ —(3lmj R ° + 2mJ *° + 2mf» } J ' (3 ' 8) 

One should note that Landau's hypothesis for quasiparticles, Hugenholtz-Van Hove theorem, the 
virial theorem and conditions of conserving approximations L '~^ ' can be exactly satisfied with the 
effective masses; the effective masses are naturally defined as the consequence of the requirement of 
self-consistency. The structure of the nonlinear mean-field approximation becomes formally equivalent 
to that of the linear a-iv Hartree approximation. 

4 Thermodynamic Consistency of the nonlinear mean-field approx- 
imation 

The Hugenholtz-Van Hove theorem and equations of motion derived by the theory of conserving 
approximations are discussed for the proof of thermodynamic consistency of the current nonlinear 



approximation . 

The energy-momentum tensor operator is given by, 

*~ = -«-/:.„ + ^j^f^y = *T + *r + ft" + t- • (4.D 

The mean-field expectation values of the baryon, scalar and vector meson operators in the ground 
state of baryons are defined by, 

d A k 



<<Br|T£» = (^(x)i^d^(x)m = -iJ j^Tr( 7 »G(kW , (4-2) 
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where the equations of motion (2.2) is used for (^/\T^ U \^/). The Roman and Greek indices in propaga- 
tors are for spin and isospin degrees of freedom; the repeated indices are summed. The momentum- 
space expressions of the eqs. (4.3) ~ (4.5) should be understood as the definition of the nonlin- 
ear mean-field approximation. The propagators, D a (k), DJ lu (k) and D .^(fe), are the full meson- 
propagators, respectively. The energy density and pressure of nonlinear a, nonlinear a-to-p mean-field 
approximations' ' ' > are all reproduced from the eqs. (4.2) ~ (4.5) by assuming corresponding 
interactions. 

The energy density, £ MFTj and pressure, Pmft, can be derived from (4.2) to (4.5), by applying 
the baryon Green's function G D (k) and the free nonlinear meson propagators with effective masses of 
mesons, D^\k) : D^ ^ v {k) and D ° a ^ v (k), to tadpole diagrams' '. One obtains the following results 
by using the relations: F MFT = (T 00 ) and p M ft = -(T tl ) (i = 1, 2, 3, i is summed) as, 

F C /^, Fm ™* 2^3, 3 ^4, 4 ml 2 9uA 4 g a 2 2 

(W 2 ) .2 / m p , fl>4 p2 , ffo-p ,2 , 9u>p T/2 \ M fi N 
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1 C f F .3; "" m o- 1 2 50-3 ,3 9cjA , 4 . ^ T ^2 , ffw4 T ,-4 . 5W ,2t/2 

p — = 3 (^f io d "^y " ~^° " it* ' it* + — y ° + ir y ° + — ° y ° 

+^Pl(f + 9 f*o + ^ + 9 -fV$) , (4.7) 

and note that mean fields of mesons are also related with baryon density as, (fro = (M — M*(k F ))/g CT , 
^0 = {9u/ m *uj')PB an d -Ro = (9p/2 mU p 2 )p3- The vacuum expectation value of the energy-momentum 
tensor must be subtracted when the finite physical energy densities are extracted from eqs. (4.3) ~ 
(4.5)! 1 ]. 

Now, thermodynamic consistency of the current approximation can be shown explicitly. By adding 
(4.6) and (4.7), one can directly calculate that the energy density, pressure and single particle energy, 
(£ MFT , Pmft, E(k)), satisfy the thermodynamic relation: 

£mft+Pmft = P B E(k F ) . (4.8) 

In order to obtain (4.8), the single particle energy, E(k) = E*(k) — TP(k F ), is used to partially integrate 
the first term of (4.6). By adding the result to the first term of (4.7), one obtains the right-hand side 
of (4.8); the rest of the terms cancel one another. The equality of chemical potential and the single 
particle energy, [i = E(k F ), is also proved by differentiating the energy density, £ M ft, with respect to 
the baryon density, p B , and using equations of motion for mesons (2.3); it results in, 

" = ^ = *<*->■ '«> 

The result (4.9) is essential for the proof of the Hugenholtz-Van Hove theorem to an approximation and 

for the definition of the density of nuclear matter saturation, where the pressure vanishes: Pmft = 0. 

The self-consistency of effective masses of mesons, m* , m* and m* , is again proved by the condition 

of conserving approximations [' 1. The functional derivative of energy density, £ M ft(4'07 ^0)-Ro> n i)> 

with respect to the baryon number distribution, rii, is given by: 
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The self-consistent condition of the theory of conserving approximations requires: — = 0, 

d(fro 

OC jut £?T ^^MF'T' rl 1 I m 

— = and — = l J . The condition independently generates meson equations of mo- 
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tion and determines self-energies of the approximation. The self-energies calculated by propagators, 
(3.5) ~ (3.7), and by the condition of conserving approximations become equivalent, only if the effec- 
tive masses of mesons are given by (3.8). The nonlinear mean-field approximation becomes thermody- 
namically consistent, relativistic, field-theoretical approximation, with the effective masses of mesons 
and coupling constants. The current functional differential approach can be applied to Hartree-Fock, 
Ring! I and other approximations, in order to obtain exact solutions which maintain the HV theorem 
and Landau's requirement of quasiparticles. 



5 Properties of nuclear matter and neutron stars 

The linear a-to Hartree approximation (LHA)[ 1, and the current nonlinear Hartree approximations, 
denoted as NHA 2 - 50 and NHA 2 - 00 which generate M max = 2.50 M© and 2.00 M , will be compared 
in this section. They are evaluated numerically by solving equations of motion iteratively with a 
given density p B ; then, the single particle energy, energy density and pressure are calculated. Since 
the Hugenholtz-Van Hove theorem is maintained, Fermi-liquid properties of nuclear matter are un- 
ambiguously determined, and properties of neutron stars are evaluated by employing the equation of 
state for neutron matter. One should note that self-consistency defined by the HV theorem, the virial 
theorem and conditions of conserving approximations are essential to calculate properties of infinite 
nuclear and neutron matter. 

The coupling constants are adjusted in order to produce saturation properties, £ / p B — M = —15.75 
MeV, at k F = 1.30 fm _1 , the symmetry energy 04 = 35.0 MeV, and simultaneously, the minimum 
value of compressibility, K, is investigated to the maximum mass, M max = 2.50 M Q and 2.00 Mq, 
of neutron stars L ' '. The calculation is repeated until required values are obtained. Contrary to 
an expectation, the adjustment of nonlinear coupling constants are not free at all because they are 
restricted by the coupled equations of motion and constraints of thermodynamic consistency; the 
values of nonlinear coupling constants are confined within a small region. Consequently, the equation 
of state and properties of nuclear matter compatible with the mass of hadronic neutron stars will be 
explicitly constrained by the density-dependent effective masses and coupling constants. The results 
should be analyzed quantitatively by nuclear experimental data in order to examine the validity of 
nonlinear mean-field models. 

The coupling constants which increase or decrease the values of both K and 04 will have the same 
tendency to the value of M max , and it depends on respective coupling constants how much they affect 
physical quantities. The a-uj mixing interaction term, g aiU , is an exception that decreases the value 
of K but increases the values of both 04 and M max , which indicates the 0-10 mixing interaction may 
simulate distinct many-body interactions at normal density. Even if one may assume the value of 
g a 3 as 10 2 ~ 10 3 , the g^-term produce small contributions to physical quantities which are readily 
recovered by changing g a ^ and g^- The properties of nuclear matter and neutron stars are modified 
mainly by changing g ff 4, g^A and the mixing term g auJ . The coupling constant, g ap N, contributes 
mainly to 04 and M max . 

The compressibility is calculated by 

* = 9 *i£ = **(!;)• <5 - 1) 

where jjl is the chemical potential and equal to the Fermi energy, fi = E(k F ), since HV theorem and 



Landau's hypothesis for quasiparticles are maintained exactly. The binding energy curves given by 
the linear a-to Hartree approximation (LHA) I J and the current nonlinear Hartree approximations, 
NHA 2,50 and NHA 2,00 , corresponding to the maximum masses of neutron stars M max = 2.50 M & and 
M max = 2.00 Mq, are compared in the Fig. 2. 

Coupling constants and Fermi-liquid properties of nuclear matter and neutron stars are listed in 
the Table I. The linear Hartree approximation (LHA) gives K = 530 MeV and M max = 3.09 Mq. The 
nonlinear Hartree approximation, NHA 2 ' 50 , shows that the lower bound of compressibility is about 
K ~ 340 MeV, whereas NHA 2 ' 00 exhibits that the lower bound of compressibility is about K ~ 180 
MeV. The results are quite distinctive that it should be investigated by studying experimental data 
of nuclear and neutron matter. 

The effective mass of nucleons, M*/M, and the effective masses of mesons, m^/rriu and m^/m^, 
of NHA 2 ' 50 and NHA 2 ' 00 in nuclear matter are shown in the Figs. 3a and 3b by comparing with the 
linear a-uj Hartree approximation (LHA). The effective masses of nucleons and mesons are slightly 
increased at saturation; the results of NHA 2 ' 50 are M* jM = 0.62, m* a /m a = 1.06 and m* /m w = 1.05 
for K = 340 MeV, a 4 = 35.0 MeV and M max = 2.50 M . The NHA 2 - 00 generates: M*/M = 0.72, 
m%/m a = 1.20 and m^/m^ = 1.17 for K = 186 MeV, a 4 = 35.0 MeV and M max = 2.00 M©. The 
effective masses of neutron matter for NHA 2 ' 50 and NHA 200 are shown in the Figs. 3c and 3d. 

Notice that the effective mass of nucleons is given by M* = M — g a (f>o and the sigma mean-field, 
<fio, increases rapidly in the LHA, resulting in the rapid decrease of M*, but the decrease of M* is 
smoothed and counterbalanced by the increase of meson masses in the nonlinear approximation. The 
density-dependences of effective constants, g* and g*, are shown in the Fig. 4. Similar to the increase of 
meson masses, the nonlinear interactions increase the values of effective coupling constants. However, 
nonlinear interactions mainly contribute to physical quantities through effective masses of baryons 
and mesons. The ctujN effective coupling constant changes only less than 3%: g^/g^ = 1.01 ~ 1.03 at 
saturation density, when the empirical data are altered from NHA 2 ' 50 to NHA 2 ' 00 . This shows that 
coupling constants are insensitive to nonlinear interactions in order to determine experimental values in 
the level of Hartree approximation. The apN effective coupling constant changes g* p /g p = 1-07 ~ 1.14 
at saturation density. The large values of g* p /g p are needed in order to reproduce the value of symmetry 
energy, a 4 = 35.0 MeV, in the Hartree approximation. However, the important contributions to the 
symmetry energy also come from the exchange interactions' ' , which will compensate the value of 
the symmetry energy and may lead to the negligible value of g p /g p at saturation density. 

The nonlinear interactions constrain effective masses of baryons and mesons, M* , m*, m* , m* 
and effective coupling constants, $*, g* Especially, the results should be compared with the rest 
masses of hadrons (M, m u , m p , m n , • • •) in order to check nuclear experimental data whether masses 
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of mesons are modified inside of nucleus or not. The quantitative study of effective masses by nuclear 
experimental data is essential to conclude the validity of nonlinear mean-field models. 
The symmetry energy is calculated by 



1 

a 4 = -Pe 



d 2 £ 



PB 



(5.2) 

P3=0 



M.. 

where p% is the difference between the proton and neutron density: p% = p p — p n = (kp — kp )/37r 2 
at a fixed baryon density, p B = p p + p n = 2k F /37r 2 . The p- meson mean-field Rq and self-energies are 
calculated by, 

g* g* 2 

-J 1 *® = ~Kp = S P« = 4~^2 P3 > ( 5 - 3 ) 

where g* p = g p + 2g up N(f)Q. The effective mass of p-meson is given by 

and it is self-consistently determined by the eq. (5.3). In symmetric nuclear matter (p$ = 0), the 
p-meson mean-field, Rq, a-p and u-p interactions vanish, and m* = m p will be restored. The effective 
mass of p-meson is similar to the bare mass: m* ~ m p at saturation density, since contributions 
to the effective mass of p-meson from nonlinear interactions are negligible in the level of mean-field 
approximation. Although the coupling constant, g P4 , is fixed from the beginning according to the 
ref. [3], one can see from the eq. (5.4) that the nonlinear self-interaction of p-meson will not contribute 
to physical qunatities at all, since the contribution to m* from the g^-term in (5.4) is negligible. The 
symmetry energy, 04, versus Fermi-momentum, k F , is shown in the Fig. 5. Nonlinear interactions 
restricted by self-consistency and properties of nuclear and neutron matter do not quite change the 
value of symmetry energy. 

It is difficult to reproduce required values of M max and a^ in the numerical calculation, since 
properties at nuclear matter saturation, coupled equations of motion and self-consistency will confine 
adjustable parameters strictly. Coupling constants should be carefully adjusted so that one can obtain 
solutions in the density range, < k F < 3.0 fm _1 , where the equation of state is important to 
determine the mass of neutron stars. The conditions of conserving approximation and consequently, 
self-consistent condition, d£/drii = E(ki), are essential to evaluate eqs. (5.1) ~ (5.2) of compressibility 
and symmetry energy. 

The equations of state for LHA, NHA 250 and NHA 200 of neutron matter are shown in the Fig. 6. 
The masses of neutron stars are calculated by Tollman-Oppenheimer-Volkoff (TOV) equation! ' ^ 
and shown in the Fig. 7, as a function of a central energy density, £ c - The data of NHA 200 indicate 
that the equation of state which produces compressibility, K > 180 MeV, is necessary in order to 
reproduce all the observed masses of hadronic neutron stars; K ~ 180 MeV is the lower bound to 
support all the hadronic neutron stars in the nonlinear a-uj-p mean-field approximation. 
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The moment of inertia, 1, of neutron stars is calculated by[ ' 1, 

T _ 8ir fR r 4 (g(r)+p(r)) e^V 2 

3 Jo (1 - 2GM(r)/r)V2 ' (5 ^ 

where G is the gravitational constant. The physical constants are set as h = c = Mq = 1, and so, 

G = GMq/c 2 = 1.476 km; the moment of inertia is given by the unit of MQ-km 2 . The metric function, 

is(r), is given by 

dv{r) _ 2GM(r) + ATrr 3 Gp(r) 

~lh~ ~ r{r - 2GM(r)) ' ^ ' 

with dv(r)/dr — ► (r — > 0), and the boundary condition is 

i/(i2) = log(l - 2GM{R)/R) , (5.7) 

at the surface of a neutron star, r = R. The masses of neutron stars and moments of inertia are 
shown in Fig. 8, Fig. 9 and Fig. 10. The approximation, NHA 250 , gives, K = 340 MeV and a^ = 35.0 
MeV, the neutron star mass, M max = 2.50 Mq, the moment of inertia, I = 368 M Q km , and the 
radius, R = 13.2 km; that is, the lower bound of compressibility is expected to be K ~ 340 MeV 
(M max > 2.50 M ). The approximation, NHA 2 - 00 , generates K = 186 MeV and a 4 = 35.0 MeV, the 
neutron star mass, M max = 2.00 Mq, I = 192 MQkm and R = 11.9 km. 

The lower bound of compressibility changes about 50% {K = 340 — ► 170 MeV), as the maximum 
mass of neutron stars is changed 20% (M max = 2.50 — > 2.00 Mq). Compared K with 04, the 
compressibility changes about 20% (K = 340 — > 270 MeV), as the symmetry energy, 04, changes 
about 8% (35.0 — > 32.5 MeV). The properties of nuclear and neutron matter are closely related 
to each other; the accumulation of data and accurate measurement of experimental values are very 
important to understand many-body systems of hadrons at nuclear density, as well as systems at 
high-energy and high-density. 

All the results are consistent with the analysis of acceptable configuration of the central energy 
and maximum masses of neutron stars^ >. The approximations LHA, NHA 2 ' 50 , and NHA 2 ' 00 produce 
distinct results for properties of nuclear and neutron stars, which may be distinguished in the course 
of the accurate examinations of nuclear experimental data and accumulating data of neutron stars. 
Since properties of nuclear matter and neutron stars are closely related, the theoretical analysis of 
self-consistency to an approximation is essential to calculate physical quantities. One should note that 
retardation corrections may break self-consistency to an approximate calculations and cause serious 
discrepancy in self-consistent calculations^ ' . The current technique of constructing a conserving 
approximation gives a method to extend an approximation self-consistently. 
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6 Conclusions and Remarks 

We have employed a nonlinear mean-field model in order to simulate complicated many-body interac- 
tions of hadrons and studied Fermi-liquid properties of nuclear matter and neutron stars. Nonlinear 
interactions renormalized self-consistently with thermodynamic consistency result in the introduction 
of effective masses and effective coupling constants of nucleons and mesons. The contributions of non- 
linear interactions to physical quantities (effective masses, M*, m*, m* , m* compressibility, K and 
symmetry energy, 04) and neutron stars (maximum mass M max , central energy density £ c , moment of 
inertia /, and radius R) are discussed and mutual relations of the physical quantities are examined. 
The nonlinear mean-field models seem to simulate properties of nuclear and neutron matter reason- 
ably well ^ ^ . However, since the effective masses of mesons and nucleons are self-consistently confined 
in order to maintain HV theorem, thermodynamic consistency' ' >, properties of infinite matter 
and coupled equations of motion for mesons, the values of nonlinear coupling constants are not free 
parameters to adjust at all. 

The nonlinear interactions are constrained with self-consistency and generate rather restricted 
values of properties of nuclear and neutron matter. The softer equation of state will decrease the 
value of compressibility, K, and simultaneously, decrease the maximum mass of neutron stars, M max 
and symmetry energy, 04 . Nonlinear interactions will connect physical quantities closely to one another 
that it is difficult to reproduce experimentally desired values of K, 04 and M max . A fine adjustment 
of nonlinear coupling constants is needed to produce the experimentally required values of K, 04 and 
M max . Therefore, it is important to examine quantitatively by nuclear experimental data, whether 
density-dependent effective masses (M*, m* , m* , m*) and coupling constants exist or not. The results 
will clarify the validity of nonlinear interactions and the relation between Fermi-liquid properties of 
nuclear and neutron stars. 

The nonlinear self-interactions of g a t± and g W 4, mixing- interaction of g ffU) terms are especially impor- 
tant to determine the binding energy, compressibility and the mass of neutron stars, but the nonlinear 
interactions yield little contribution to symmetry energy in the level of Hartree approximation. The 
density-dependence of effective coupling constants, 5* and g*, are not so important to reproduce the 
symmetry energy 04 and M max of neutron stars. The existence of solutions to coupled equations of mo- 
tion for mesons and convergences of properties of nuclear matter in the density range, k F = 1.0 ~ 3.0 
fm _1 , should be carefully checked in order to calculate properties of neutron stars. The values of non- 
linear coupling constants might be studied as renormalized values of many-body interactions, such as 
Hartree-Fock, Ring, Ladder and other effective approximations with vertex interactions; the effective 
masses and coupling constants should be investigated further in terms of many-body interactions. 

A mixed-phase of quark-hadron star might be expected. A hadron-quark matter mixed-phase 
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equation of state which is defined by MIT-bag model, for example, could produce the maximum 
observed mass: M max = 2.50 M@[ '. However, even if the experimental value of compressibility and 
the maximum mass of neutron stars are expected to be K ~ 200 MeV^ 1 and M max ~ 2.00 M Q , 
the nonlinear a-u-p mean-field approximation indicates that it is possible to generate neutron stars 
which might be pure hadron-stars. In order to distinguish the inner structure of neutron stars from 
the observed experimental data, it is important to analyze effective masses of baryons and mesons 
from nuclear experimental data as precisely as possible. 

The current approach by employing the theory of conserving approximations is essential to check 
self-consistency and calculate equations of state, Fermi-liquid properties and Landau parameters^ > 
of nuclear matter. It provides us with a method to construct thermodynamically consistent approxi- 
mations. The nonlinear mean-field approximation suggests that nonlinear interactions of mesons and 
baryons be perceived as effective masses and effective coupling constants. These renormalized con- 
stants are essential and naturally required for self-consistency to an approximation and might be used 
to examine the validity of nonlinear models of hadrons. The accumulation of data and the analysis of 
neutron stars should be investigated further in order to elucidate the validity of nonlinear mean-field 
models and theory of nuclear models. 
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Table I. Properties of nuclear matter and neutron stars 

The results of the linear a-co Hartree approximation (LHA)I 1, and the current nonlinear Hartree 
approximations, NHA 2 - 50 and NHA 2 - 00 for M max = 2.50 M Q and 2.00 M , are shown respectively. 
M* is the effective mass of nucleon; to*, to* and to* are effective masses of a, to and p mesons, 
respectively. K is compressibility and 04 is symmetry energy. The masses are fixed as M = 939 MeV, 
m a = 550 MeV, to w = 783 MeV and m p = 770 MeV^ 1 . M max is the maximum mass in the solar mass 
unit (Mq) and Sc (10 15 g/cm 3 ) is the central energy density; / is the inertial mass (Mq ■ km 2 ) and 
R (km) is the radius of a neutron star. The coupling constants are adjusted to reproduce the values: 
£ / p B — M = —15.75 MeV, at k F = 1.30 fm _1 , and 04 = 35.0 MeV; the maximum masses of neutron 
stars, M max = 2.50 Mq and M max = 2.00 Mq 1 >i >. The values of vertex coupling constants are, 
9aujN = 0.0050 and g apN = 0.0080 for NHA 2 - 50 ; g awN = 0.0120 and g apN = 0.0215 for NHA 200 . The 
coupling constants, g p and g p 4, are fixed' 1 J - ^. 





9a 


9u> 


9p 


9*3 (MeV) 


9oA 9wA 


9pi 


9cru> 


9ap 


9uip 


LHA 


11.071 


13.795 


8.077 


0.0 


0.0 0.0 


0.00 


0.0 


0.0 


0.0 


NHA 2 - 50 


10.647 


12.863 


8.077 


100 


280 280 


4.00 


70 


2.0 


2.0 


NHA 2 - 00 


10.263 


11.737 


8.077 


480 


1700 1700 


4.00 


440 


21.0 


21.0 


M*/M m%/m a 


mllm w 


K (MeV) 


a 4 (MeV) 


Mmax 


Sc 


I 


R (km) 


LHA 


0.54 


1.00 


1.00 


530 


35.0 


3.09 


1.25 


717 


14.5 


NHA 2 - 50 


0.62 


1.06 


1.05 


340 


35.0 


2.50 


1.53 


368 


13.2 


NHA 2 - 00 


0.72 


1.20 


1.17 


186 


35.0 


2.00 


1.94 


192 


11.9 
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Figure Captions 

Fig. la. The scalar self-energy, £ s , is given by renormalized effective masses and nonlinear coupling 
terms. The vertices are — ifiWivVcn for era; iV- vertex and —ig a pNRoT3'J for apiV-vertex. 

Fig. lb. The omega and rho meson sources to generate £ . The vertex factors are —ig^l^ f° r uN 
and — i(<7*/2)r37 M for pN couplings. 

Fig. 2. The binding energy curves of symmetric nuclear matter. The linear a-to Hartree approxima- 
tion (LHA), and the nonlinear Hartree approximations, NHA 2 ' 50 and NHA 2 ' 00 which produce 
the maximum masses of neutron stars: M max = 2.50 M and M max = 2.00 Mq, are shown 
respectively (see, the table I). 

Fig. 3a. The effective nucleon mass M*/M (solid-line), the effective scalar meson mass, m* s /m s 
(dotted- line), and the effective vector meson mass, m^/rn^ (dot-dashed line), for the binding 
energy of nuclear matter, NHA 2 - 50 (K = 340 MeV, M max = 2.50 M & ). 

Fig. 3b. The effective masses of nucleons and mesons for nuclear matter, NHA 200 {K = 186 MeV, 
M max = 2.00 M ). 

Fig. 3c. The effective masses of nucleons and mesons for neutron matter, NHA 250 (K = 340 MeV, 
M max = 2.50 Mq). The degeneracy factor is £ = 2, with the same coupling constants as those 
of Fig. 3a (see, the Table I). 

Fig. 3d. The effective masses of nucleons and mesons for neutron matter, NHA 200 (K = 186 MeV, 
M max = 2.00 M & ). The degeneracy factor is C = 2, with the same coupling constants as those 
of Fig. 3b (see, the Table I). 

Fig. 4. The effective coupling constants, g^/9w (solid line) and g* p /g p (dashed line). In the solid and 
dashed lines, the lower lines are for NHA 250 and the upper lines for NHA 2,00 . 

Fig. 5. The symmetry energies, LHA, NHA 2 ' 50 and NHA 2 ' 00 , versus Fermi-momentum, k F , for Table I 
are shown. 

Fig. 6. The equations of state for neutron matter, p = S is the causal limit of the equation of state. 

Fig. 7. Neutron star masses as functions of the central energy density, £ c . The lines are LHA {M max = 
3.09 M Q , £ c = 1.25 x 10 15 g/cm 3 ), NHA 2 - 50 {M max = 2.50 M , £ c = 1.53 x 10 15 g/cm 3 ), and 
NHA 2 - 00 {M max = 2.00 M , £ c = 1.94 x 10 15 g/cm 3 ), respectively. 

Fig. 8. The masses and radii of neutron stars. 
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Fig. 9. The moments of inertia and radii of neutron stars. 
Fig. 10. The moments of inertia and masses of neutron stars. 
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Fig. 10 
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